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Abstract 

We examine the corecoUapse times of isolated, two- 
mass-component star clusters using Fokker-Planck 
models. With initial condition of Plummer models, 
we find that the corecoUapse times of clusters with 
M1/M2 ^ 1 are well correlated with 
{Ni/N2y/^{mi/m2)^trh, where M1/M2 and 7711/7712 
are the light to heavy component total and individ- 
ual mass ratios, respectively, N1/N2 is the number 
ratio, and trt is the initial half-mass relaxation time 
scale. We also find two-component cluster parameters 
that best match multi-component (thus more realis- 
tic) clusters with power-law mass functions. 
Subject headings: celestial mechanics, stellar dynamics — 
globular clusters : general 

1. INTRODUCTION 

The course of dynamical evolution of pre- and post- 
collapse globular clusters is determined by many fac- 
tors such as initial mass function, nature and effi- 
ciency of energy generation mechanisms, tidal cut-off, 
anisotropy of velocity distribution, initial population 
of binaries, and stellar evolution. There have been 
many efforts in developing more and more complex 
cluster models including such factors, making analy- 
sis and interpretation rather difficult. To study the 
dynamical evolution of globular clusters more real- 
istically, among others, Chernoff & Weinberg (1990) 
included the effects of stellar evolution, Lee, Falhman, 
& Richer (1991) used multi-component models, Taka- 
hashi (1995) includeded an anisotropic velocity distri- 
bution. 
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However, studying simpler models could be more 
instructive in identifying important physical processes 
governing the evolution. Kim, Lee, & Goodman 
(1997; hereafter KLG) studied on the postcollapse 
evolution of cluster variables and the gravitational os- 
cillation using two-component Fokker-Planck models. 
In this paper, as a supplementary study to KLG, we 
present a fitting formula for the corecoUapse times 
of two-component models and compare the results of 
two- and multi-component models. As in KLG, here 
both tidal-capture binary heating and tree-body bi- 
nary heating are included, and clusters are assumed 
to be isotropic and isolated (no tidal cutoff). For 
the methods that we are using here and the benefits 
of studying simpler models (two-component models), 
readers are referred to KLG and references therein. 

CorecoUapse times of two-component clusters were 
presented by Inagaki & Wiyanto (1984), Inagaki (1985), 
and Lee (1995) among others. These papers calcu- 
lated corecoUapse times as a function of M2/M1, the 
ratio of total masses of heavy component to light com- 
ponent, and found that the ratio of corecoUapse time 
to initial half-mass relaxation time, tcc/trh, has a min- 
imum value at M2/M1 ^ 0.1. However, the paramter 
M2/M1 may be divided into 7772/7)11, the individual 
mass ratio, and N2/N1, the number ratio. In the 
present paper, we calculate the corecoUapse times of 
two-component models as a fuction of more complete 
two-copmonent cluster parameters, and find a fitting 
formula between them. However, the clusters stud- 
ied here are restricted to those with Mi 3> M2 as in 
KLG. 

On the other hand, it would be helpful in inter- 
preting the results of two-component models if the 
similarities and discrepancies between the results of 
two- and multi-component models are well known. In 
the present paper, we also compare two-component 
models to 11-component models, and thus provide a 
way to extrapolate to more realistic cluster the results 
of two-component clusters such as those in KLG. 

2. CORE-COLLAPSE TIMES 

To calculate the corecoUapse times of two-component 
star clusters, we have performed total 11 runs of direct 
numerical integration of the orbit-averaged Fokker- 
Planck equation with a local approximation. The 
code used here is descended from Cohn (1980). Pa- 
rameters of our two-component runs are shown in Ta- 
ble n^. This set of parameters has been chosen such 



that it provides all possible combinations of parame- 
ters M, N, 1712/ nil, and N1/N2, where M is the clus- 
ter mass and N is the total number. Note that in all 
our runs, the total mass of heavy component, M2, is 
negligible compared to the total mass of light compo- 
nent. Ml, and thus mi ~ M/Ni. The initial density 
and velocity profiles are given by Plummer models 
with Vci/vc2 = 1 and Pci/pc2 — M1/M2, where is 
the three-dimensional core velocity dispersion, and pc 
is the core density. 

CorecoUapse times of our runs are shown in Table 
in units of 10^° yr and trh- We empirically found 
that tec can be fitted by the following formula: 
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where N5 = iV/10^ M5 = M/lO^Mg, and rn is the 
initial half-mass radius. Each tec value is plotted over 
the righthand side of the above equation in Figure 
which shows a good X-Y correlation. Equation 
) is to be compared with the standard half-mass 
relaxation time scale, 
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where Vm is the root-mean-square three-dimensional 
velocity of the whole cluster and (w^)i>=t>„ is the aver- 
age change of in parallel component to initial Vm 
per unit time. 

Isolated single-mass clusters with initial condition 
of Plummer models collapse at 15.4 trh (Cohn 1980), 
where trh is the half-mass relaxation time scale and 
does not vary much until the corecoUapse takes place. 
However, the ratios of the time required for corec- 
oUapse tec to trh and core relaxation time scale trc 
strongly depend on the density and velocity profiles. 
Quinlan (1996) found that for single-mass clusters, tec 
varies much less when expressed in units of trc divided 
by a dimensionless measure of the temperature gradi- 
ent in the core. Although in single-mass clusters the 
velocity profile (as well as other physical parameters) 
evolves by the two-body relaxation, two- or multi- 
component clusters have another driving force: the 
equipartition. 

Both mass segregation and equipartition are en- 
volved in determination of the time to corecoUapse, 
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Fig. 1. — CorecoUapse times of runs in Group A. Nr, = 
and M5 = Mq. 

and this complexity makes the theoretical interpreta- 
tion of the above correlation between tec and cluster 
parameters quite difficult. Here we suggest the follow- 
ing analysis as one way to explain this correlation. 

The actual duration of corecoUapse is very small 
compared to the time to corecoUapse from the be- 
gining of cluster's evolution. Instead, clusters spend 
most of their precollapse phases under mass segrega- 
tion process and approach to the onset of homolo- 
gous phase of corecoUapse. If a considerable degree of 
equipartition is accomplished in the precollapse phase 
as in all of our two-component models, the time to the 
onset of corecoUapse will be determined by how fast 
the light component gains the energy from the heavy 
component via equipartition. Thus one may define 
the precollapse time scale of two-component clusters 
tr2 as following: 



(3) 



where (t'p)2 is the velocity dispersion change of the 
light component via interactions with heavy compo- 
nent. Using the standard expression for the average 
velocity dispersion change per unit time, one has 



^1)2 oc 



G'^m2Pm2 



Vm2 



(4) 



2 



Table 1 

Parameters and Core-Collapse Times of Two-Component Models 



Run 


1712 

mi 


Ni 

N2 


M 
(Mo) 


iV 


TO2 

(Me) 


^CC 

(lOiOyr) 


^cc / ^rh 


baab 


2 


100 


iqs 


141457 


1.4 


1.27 


12.42 


caab 


3 


100 


iqs 


210125 


1.4 


0.93 


6.34 


faab 


4 


100 


iqs 


277473 


1.4 


0.61 


3.23 


cdab 


3 


30 


iqs 


201299 


1.4 


0.56 


3.97 


cbab 


3 


300 


iqs 


212871 


1.4 


1.62 


10.92 


caabl 


3 


100 


10^ 


70042 


3x1.4 


0.35 


6.47 


caab2 


3 


100 


10^ 


630374 


ixl.4 


2.52 


6.28 


baabS 


2 


100 


iqs 


212185 


|xl.4 
|xl.4 


1.80 


12.17 


faabS 


4 


100 


iqs 


208104 


0.47 


3.23 


caeb 


3 


100 


3 X W 


63037 


1.4 


0.57 


6.34 


cabb 


3 


100 


3 X 10^ 


630374 


1.4 


1.49 


6.43 



Note. — The initial half-mass radii rh of these runs are all 5 pc. 



equation, we find that the degrees of equipartition 
in the precollapse phases of our two-component runs 
do not directly correlate with the minimum values of 
equation (|^). In fact, exact equipartition is usually 
not accomplished even when the value of equation 
is less than unity, because as mass segregation 
of heavy component progresses, interactions between 
heavy and light components occur less. Thus equa- 
tion (H) should be regarded as a degree of tendancy to 
equipartition and it is this tendancy that tr2 requires 
in its definition. 

While Quinlan (1996) introduced a temperature 
gradient in the core in a derivative form to explain a 
huge variation in tec for clusters with different initial 
profiles, here we introduced both density and velocity 
gradients of heavy component naturally into the time 
scale by considering global equipartition. 

3. COMPARISON WITH MULTI- 
COMPONENT CLUSTERS 

Clusters have continuous mass functions. How- 
ever, mass functions are usually realized with dis- 
crete mass components in numerical calculations. Sci- 
entists found that 10 to 20 components are enough 



where the heavy component mean density pm2 is pro- 
portional to M2/r^ and the Coulomb logarithm has 
been omitted. It is also assumed that v^^i ~ fm2- 
Spitzer (1969, 1987) showed that for a two-component 
cluster of polytropic index n between 3 and 5 with 
Ml ^ M2 and a Maxwellian velocity distribution in 
a parabolic potential well, the minimum degree of the 
global equipartition is a function of cluster's parame- 
ters such that 
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With the minimum value of equation (^ and as- 
sumptions that Ml ^ M2 and ~ GM/rh, equa- 
tion (^ now becomes 

5/3 

tr2 ^ i'-^] { — ] NM-'/\'J^ 



Ni' 


2/3 


mi 


N2, 


) 


Km2 


Ni' 


2/3 


( nil 


N2, 


) 


\m2 



- l^j Kit)"'-'- 

The above is in the same form as equation (|l|) with 
only small discrepancies in the exponents. Although 
equation (||) has been used for derivation of the above 
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to represent continuous mass functions for Fokker- 
Planck models, and such numerical representation 
for a given mass function is quite straightforward for 
these multi-component clusters: there is only a ques- 
tion of choice of each component's mass bin and a 
representative value. However, when the number of 
components is reduced to 2 for the sake of analytical 
simplification, such choice is not so simple because 
dynamically important mass and corresponding num- 
ber of stars may be different from simple mean mass 
and total number of a certain mass range. Therefore 
two-component cluster parameters (such as 7712/7711, 
N1/N2, and N) that well represent a continuous mass 
function should be numerically found through com- 
parisons of the evolution of two- and multi-component 
clusters. 

In this section we will compare our multi-component 
models with the two-component models in KLG vary- 
ing M and the mass function of the multi-component 
models. Cluster parameters of our multi-component 
models are given in Table The initial density and 
velocity profiles are given by Plummer models. The 
initial half-mass radii of all multi- and two-component 
models are 5 pc. The number of component is 11 and 
we adopt a power-law mass function: 



N{m)dm oc 777 '^^"'"■'dTn, 



(7) 



where x is the mass spectral index and the Salpeter 
mass function has x = 1.35. For a bin i with bound- 
aries rriia and mn,, the total mass in the bin is ob- 
tained by 



niN{m)dm. 



(8) 



Then the number of stars in the bin is Ni — Mi/rrii, 
where 77ii is the representative mass of each bin. The 
main-sequence star mass range was selected to be 
0.08 M0 - 0.8 Mq. Following Sigurdsson & Phinney 
(1995), the stars of initial mass nii between 0.8 Mq 
and 4.7 Mq were assumed to have evolved to white 
dwarfs of mass 0.58 -I- 0.22 x (777,ms — 1-0) Mq, where 
Toms is the main-sequence mass, while stars of 771ms 
between 4.7 Mq and 8.0 M0 were assumed to disrupt 
completely. The stars heavier than 8.0 Mq but lighter 
than 15.0 Mq were assumed to become neutron stars 
of mass 1.4 Mq. Neutron stars are born with a kick 
velocity due to an asymmetric explosion, and they are 
ejected from the cluster if the kick velocity is greater 
than the escape velocity of the cluster. However, we 
assumed that all neutron stars remain in the cluster, 
because the retention rate of neutron stars are not 



well known and the precise realization of real clusters 
is not our goal in this study. The mass range, repre- 
sentative mass, and number of stars of each compo- 
nent is shown in Table ^. Our multi-component mod- 
els include both three-body binary heating and tidal- 
capture binary heating, but we find that the postcol- 
lapse phases of all our runs are driven by three-body 
binary heating. 

We find a two-component model which best de- 
scribes a given multi-component model by comparing 
the values of cluster variables pc, Vc, a-t t = 10^^ yr, 
and tec- An epoch of 10^^ yr has been selected as in 
KLG because by that time, our runs have reached 
self-similar expansion phase. With two-component 
models, KLG found the following numerical values: 



Pc 



4.5 X 10^ Mq 



-10/3 



;vflO/3,-2.0. 



(9a) 



V, ^ 3.8km/s (J^^ '^\l/Hl}Jh^°-^^; (9b) 



rh 



I "12 

0.042 pc — 
Vnii 



7/6 



^-4/3j^^l/3^a65. (9^) 



(9d) 



where = iV/10^ Afg = M/IO^Mq, and in ee 
t/lO^^yr. On the other hand, the numerical values 
from our multi-component models are given in Ta- 
ble ^ There are four two-component parameters to 
be determined for a given multi-component model, 
777.2/7(71, N1/N2, N, and M. However, since clus- 
ter variables at a certain epoch in the postcollapse 
phase are independent of N1/N2 as in equation (^, 
N1/N2 has to be determined from the corecoUapse 
time, equation (|l|). Then the rest three parameters, 
777.2/777.1, N, and M, may be determined from equa- 
tion (^. This method will be called Method A, and 
parameters obtained in this way are given in Table ^ 
With Method A, our multi-component model B2 
is best described by a two-component model with 
m2/mi = 2.3, N1/N2 = 29, N = 1.4 x 10^ and 
M = 0.97 X 10^ Mq. Note that with these parameters, 
77^2 ~ 1.6 Mq, which is little higher than the mass 
of the heaviest component of our multi-component 
models, 1.4 Mq. Neutron stars play an important 
role in dynamical evolution of globular clusters: a 
considerable fraction of dynamical binary formation 
(as apposed to primordial binaries) envolves neutron 
stars. For this reason, in two-component clusters, the 
heavy component is often targeted for neutron stars 
and the light component is for main-sequence stars. 
Thus in finding the best matching two-component pa- 
rameters, it could be more meaningful if 777,2 is set to 
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Table 2 

Parameters and Results of Multi-Component Models 



Values at t = 10^^ yr 



Run 


X 


M 


N 




Pc 


Vc 












(lo^V) 


(Mqpc-3) 


(kms-i) 


(pc) 


A2 


1.00 


10^ 


402857 


0.528 


7.63 X lO'' 


2.79 


32.8 


B2 


1.35 


10^ 


509201 


0.584 


8.13 X 10* 


2.75 


27.7 


C2 


1.50 


10^ 


552232 


0.621 


8.19 X 10" 


2.73 


26.2 


Bl 


1.35 


3 X ID" 


152760 


0.356 


1.71 X 10^ 


1.30 


38.5 


B3 


1.35 


3 X 10^ 


1527603 


0.931 


2.58 X 10^ 


5.44 


20.7 



Note. — The initial half-mass radii th of these runs are all 5 pc. 



Table 3 

Mass Spectra of Multi- Component Models 



Bin rrii Mass Range a; = 1.00 a; = 1.35 x = 1.50 





(Mo) 


(Mo) 






rrii ■ Ni 




Ni 


ni l ■ Ni 






m.i ■ Ni 


1 


0.1 


0.08 - 


0.15 


1, 


.00000 


1.00000 


1, 


.00000 


1. 


.00000 


1, 


,00000 


1. 


,00000 


2 


0.2 


0.15- 


0.25 


0, 


.40631 


0.81262 


0. 


.33263 


0, 


,66526 


0, 


,30517 


0, 


,61034 


3 


0.3 


0.25 - 


0.35 


0, 


.17842 


0.53526 


0, 


.12584 


0, 


.37752 


0, 


,10826 


0, 


.32478 


4 


0.4 


0.35 - 


0.45 


0, 


.09995 


0.39980 


0, 


.06359 


0, 


,25436 


0, 


,05223 


0, 


,20892 


5 


0.5 


0.45- 


0.55 


0, 


.06385 


0.31925 


0, 


.03753 


0, 


,18765 


0, 


,02985 


0, 


,14925 


6 


0.6 


0.55 - 


0.65 


0, 


.11712 


0.70272 


0, 


.05766 


0, 


.34596 


0, 


,04264 


0, 


.25584 


7 


0.7 


0.65 - 


0.80 


0, 


.09010 


0.63070 


0, 


.04104 


0, 


,28728 


0, 


,02945 


0, 


,20615 


8 


0.9 


0.80- 


1.0 


0, 


.02428 


0.21852 


0. 


.00822 


0, 


,07398 


0, 


,00516 


0, 


,04644 


9 


1.1 


1.0 - 


1.2 


0, 


.01286 


0.14146 


0, 


.00389 


0, 


.04279 


0, 


,00232 


0, 


.02552 


10 


1.3 


1.2 - 


1.4 


0, 


.00774 


0.10062 


0, 


.00215 


0, 


,02795 


0, 


,00124 


0, 


,01612 


11 


1.4 


1.4 - 


1.4 


0, 


,00913 


0.12782 


0. 


.00182 


0, 


,02548 


0, 


,00091 


0, 


,01274 



Note. — Ni and m, • Ni are normalized with bin 1 values. 
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Table 4 

Best Matching Two-Component Cluster Parameters 



Method A Method B 



Run 


m2 
mi 


Ni 

N2 


N5 


Ms 


m2 
mi 


Ni 
N2 


N5 


Ms 


^2 


A2 


1.78 


13.1 


1.05 


0.90 


1.70 


10.9 


1.00 


0.82 


8400 


B2 


2.33 


28.5 


1.40 


0.97 


2.17 


21.3 


1.30 


0.84 


5800 


C2 


2.56 


39.4 


1.54 


0.99 


2.36 


28.4 


1.41 


0.84 


4800 


Bl 


2.74 


53.9 


0.51 


0.35 


2.37 


30.0 


0.44 


0.26 


1400 


B3 


2.05 


18.1 


3.47 


2.49 


2.00 


16.3 


3.38 


2.36 


20000 



Note. — N5 = N/IO^ and M5 = A// 10^ Mq. N2 has been approximately 
calculated by N/{Ni/N2 + 1) and has only two significant digits. 



1.4M0. With a restriction of m2 = 1.4Mo, now 
the number of variables in equation required for 
determination of two-component cluster parameters 
is reduced to two. Since pc and are two cluster 
variables that represent the status of the core and en- 
velope, respectively, we use these variables along with 
7712 = 1.4 Mq and equation (Q) for our second method 
(Method B) to find the best matching two-component 
model (see Table H). 

With Method B, model B2 is now best described by 
a two-component model with 7712/7711 = 2.2, N1/N2 = 
21, iV = 1.3 X 10^ and M = 0.84 x lOHf©- Note 
that with these parameters, A^2 = 4890 and this value 
is about the same with the number of stars in the 
heaviest four bins (bins only for degenerate stars) of 
model B2 (X]i=8 11 ~ 5176). This may imply that 
the epoch of corecollapse is mainly determined by the 
number of stars above the turnoff mass. This fact also 
holds for other runs with different M and x. 

For clusters with N (x M (as for our multi- 
component clusters Bl, B2, and B3), equation (^) 
may be written as pc oc M^^^'^, Vc oc M^/"^, and 
rh oc From runs Bl, B2, and B3 in Table |, 

Pc Vc, and rfi are found to be proportional to M'^'^^, 
M°-622^ ^j^j Af-"-269, respectively. The absolute val- 
ues of these exponents are little smaller than equa- 
tion (^. However, since the discrepancies are not 
so significant, we conclude that the evolution aspects 
of the postcoUapse multi-component clusters are still 



well predictable from the numerical and analytical re- 
sults of two-component clusters. On the other hand, 
cluster variable values at t — 10^^ yr show relatively 
small X dependence. 

The results from Method B in Table ^ indicate 
that multi-component clusters may be described by 
two-component clusters with masses 15 to 20 % less 
and with r?7i near the turnoff mass. Of course, this 
TTii is dependent on x such that clusters with steeper 
mass function are matched by two-component clus- 
ters with smaller ttii. However, interestingly, lighter 
multi-component clusters with the same x also re- 
quire smaller 7771. This comes from the fact that 
N (X M holds for runs Bl, B2, and B3, while not 
for their matching two-component clusters by Method 
B (A^ oc M°-^^). For Method A, the best matching 
two-component clusters of runs Bl, B2, and B3 show 
A^ oc which results in nearly the same mi. 

Thus we conclude that the above difference in 7771 val- 
ues by Method B for clusters with the same x stems 
from the restriction, 7772 = 1.4 Mq. 

The evolution of the two-component model with 
the above parameters is plotted in Figure || as well as 
that of model B2. Cluster variables pc, Vc and r/j of 
two runs well coincide. Only the corecollapse times 
show a small discrepancy. This is partly because of 
the dispersion of tec from the fitting formula, equation 
(|l|), and partly because of the small N1/N2 value: 
equation (|l|) is to be used for clusters with Mi ^ M2 ■ 



6 




10 

log(t/yr) 

Fig. 2. — Comparison of the evolution of multi- 
component model (run B2; thick lines) and best- 
matching two-component model by Method B 
(m2/mi = 2.2, iVi/Ar2 = 21, N = 1.3 x 10^ and 
M = 0.84 X 10^ Mq; thin lines). The units of Pc, Vc, 
and rh are MqPc"^, kms~^, and pc, respectively. 



4. SUMMARY 

Wc have investigated the evolution of isolated two- 
component clusters with initial condition of Plummer 
models. The corecoUapse time tec showed a good cor- 
relation with a parameter (7Vi/A^2)'^^^(™i/™2)^ir/i- 
To explain this correlation, a new time scale for 
the precollapse evolution of two-component clusters, 
tr2 = {Ni/N2Y^'^{mi/m2)^^^trh have been intro- 
duced using Spitzer's (1969, 1987) global equiparti- 
tion analysis. 

We also found two-component clusters which best 
match with our multi-component clusters with power- 
law mass functions. For example, the evolution of 
11-componcnt cluster with a Salpeter mass function 
and M = 10^ Mq was well described by a two- 
component cluster with TO2/mi = 2.2, N1/N2 = 21, 
AT = 1.3 X 10^ and M = 0.84 x 10'"^ Mq. Fur- 
thermore, it has been found that the best matching 
two-component cluster has N2 very close to the num- 
ber of stars heavier than turnoff mass of the multi- 
component cluster. 
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